Constrained fitting of three-point functions 
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We determine matrix elements for B — » D semileptonic decay. The use of the constrained fitting method and 
multiple smearings for both two- and three-point correlators allows an improved calculation of the form factors. 



1. Constrained curve fitting 

The aim of this contribution is to illustrate 
the method of constrained curve fitting by ap- 
plication to three-point functions describing the 
semileptonic decay of heavy-light mesons. For a 
discussion of the method sec Jjj. The constraint 
method offers the chance to include knowledge 
about the fit parameters and generally results in 
very stable fits. This allows for much more ambi- 
tious fits than previous methods and the results 
are less prone to systematic errors from fitting. 

We apply the method to the subset (8 = 5.7, 
a tuq — 2, 4, 8) of the data analysed in || , where 
NRQCD is used for the heavy quark. The lattice 
size is 12 3 x 24. The smearings used are a Gaus- 
sian with a radius of a tq = 2.0 for the heavy and 
local for the light quark (referred to as smear- 
ing 1) and ar q /Q = 3.0 for both (smearing '2). 
All combinations of these smearings at source and 
sink were used. It has been demonstrated that 
these are not the optimal smearings J3|, in par- 
ticular for nonvanishing momenta, which leaves 
room for improvement of the results presented 
here. 



2. Form factors 

We aim at extracting the form factors h + {u>) 
from the matrix elements 



(B\V \D) = F^q 2 ) ( Pb +Pd) 

= y/m B m D h + (uj) (v B + v D )o 

for the elastic case p B = p 2 D with degenerate 
masses ms = TOrj . Vo is the leading order heavy- 
heavy current for NRQCD {§. 

First we extract the form factors for each value 
of q 2 . We fit the two-point functions parametrised 
by 

n 

Gf =Y,AiAiexp(-E k (t-t )) 

k=l 

simultaneously with the three-point functions 
using the ansatz 

= AiJuAie-^-^e-nP-Q 

k,l=l 

with current matrix elements Jki, where i,j la- 
bel the smearings and k, I the energy levels. We 
set Jki = Jik, arid Jn is h + (u>). The distance 
between quark source and sink is 12. Instead of 
the energies we in fact used as parameters the 
logs of energy differences as in 0. Our standard 
priors are 0.7(*/2) for the ground state energy^ 
and 0.6(*/2) for the energy differences between 

*i.e. our prior favours In 0.7 ± In 2 for the fit parameter 
mf? E . s . 



2 



the excitations. The first is known from previ- 
ous fits; the excitation energy priors are appro- 
priate to typical radial excitations at this lattice 
spacing. For the amplitudes we took as priors 
1.0(±1.5) and 0.5(±1.5) for the ground states, 
-0.4(±2) for A\ and 0.1 (±2) for the rest. These 
are based on experience from previous fits but 
allow for plenty of variation (the typical errors 
on the amplitudes from the fits are 0.01 . . . 0.05). 
The priors for the Ju are 1(±0.5) for Jkk and 
0.1(±0.5) for J kl (k^l). 
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very small (as the examples in |l|), the priors 
contribute only marginally. Quoted errors are 
extracted from the curvature matrix. Figure |l| 
shows a typical fit result for the three-point func- 
tions. To shed light on the effect of the priors we 
repeated the calculation with priors away from 
the expectations and allowing for smaller devia- 
tions. We find that the data strongly constrain 
the lowest lying current matrix elements Ju and 
J12, but have little information about the higher 
ones. This is expected because off-diagonal ma- 
trix elements should be, and are, small, because 
they represent the overlap (as q 2 — > 0) of orthog- 
onal states. 

The dependence of the lowest currents on the 
number of exponentials in the fit is demonstrated 
in Fig. |^. Stability sets in for six exponentials. 
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Figure 1. Three-point functions and 6-exp fit re- 
sults (p 2 = 1, q 2 = 2, mq = 2); the gray boxes are 
(displaced) enlargements of the functions under- 
neath 



Figure 2. The current matrix elements for the 
lowest lying states vs. the number of exponentials 
used in the fit (p 2 — 1, q 2 = 2, ttlq = 2) 



3. Kinetic masses 



The minimisation was done using the 
Levenberg-Marquardt algorithm. The value of 
Xaug f° r tne reliable fits (i.e. using a sufficient 
number of exponentials) is generally about 1. 
The contribution of the two-point functions is 



To convert these form factors into h + (uj) one 
needs the kinetic masses. They are calculated 
from the energy splittings AE between the p 2 = 1 
and p 2 = states 

m=(p 2 ~(AE) 2 )/(2AE). 
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Our method permits us to extract them di- 
rectly from simultaneous fits of the respective 
smeared-smeared two-point functions. This is an 
example where the usage of the classical fitting 
method could be misleading if one considers the 
two-point functions separately, because the two- 
point functions seem to plateau much earlier than 
their ratio. 

The values for the kinetic masses for the three 
heavy quark masses investigated are listed in the 
following table: 



m Q 


AE 




2 


0.058(1) 


2.9(2) 


4 


0.029(3) 


4.7(5) 


8 


0.019(3) 


7.3 (1.3) 



This much simpler fit shows stable results al- 
ready using three exponentials (see Fig. |^) . 

In the future we plan to fit simultaneously all 
the two- and three-point functions at different 
momenta and the kinetic mass will then come di- 
rectly from the same fits that give the current 
matrix elements. 
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Figure 3. Kinetic mass as function of the number 
of exponentials in the fit; the horizontal lines arc 
the perturbative results from 0] 



4. Isgur-Wise function 

Now we can plot the form factors as function 
of u> — v ■ v' , the Isgur-Wise function (Fig. ||). 
The errors are half of those from using only one 
smearing. Each dataset extrapolates to the point 
(1,1) as it has to for NRQCD. There seems to 
be a tendency of increasing slope with the heavy 
quark mass which has to be explored further. 
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Figure 4. Form factors as functions of uj 
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